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Abstract 


The  adaptive  local  polynomial  Fourier  transform  is  employed  for  the  improvement 
of  the  ISAR  images  in  complex  reflector  geometry  cases,  as  well  as  in  cases  of 
fast  maneuvering  targets.  It  has  been  shown  that  this  simple  technique  can  produce 
significantly  improved  results  with  a  relatively  modest  calculation  burden.  Two 
forms  of  the  adaptive  LPFT  are  proposed.  The  adaptive  parameter  in  the  first  form 
is  calculated  for  each  radar  chirp.  An  additional  refinement  is  performed  by  using 
the  information  from  the  adjacent  chirps.  The  second  technique  is  based  on  the 
determination  of  the  adaptive  parameter  for  different  parts  of  the  radar  image.  The 
numerical  analysis  demonstrates  the  accuracy  of  the  proposed  techniques.  It  is 
important  to  note  that  the  proposed  technique  does  not  assume  any  particular  model 
of  radar  target  motion.  It  can  be  applied  for  any  realistic  motion  of  targets. 


Resume 


La  transformee  de  Fourier  polynomiale  locale  (TFPL)  adaptative  est  ici  utilisee 
pour  ameliorer  les  images  ISAR  dans  le  cas  de  reflecteurs  a  geometrie  complexe  et 
de  cibles  effectuant  des  manoeuvres  rapides.  II  a  ete  demonte  que  cette  technique 
simple  ameliore  sensiblement  les  resultats  au  prix  d’un  effort  de  calcul  relative- 
ment  modeste.  Deux  formes  de  la  TFPL  adaptative  sont  proposees.  Le  parametre 
adaptatif  de  la  premiere  forme  est  calcule  pour  chaque  impulsion  comprimee  radar. 
Le  resultat  est  raffine  a  l’aide  de  V information  provenant  des  impulsions  compri- 
mees  adjacentes.  La  deuxieme  technique  se  fonde  sur  la  determination  du  parametre 
adaptatif  pour  differentes  parties  de  l’image  radar.  L’analyse  numerique  montre 
l’exactitude  des  techniques  proposees.  Notons  que  la  technique  proposee  ne  sup¬ 
pose  aucun  modele  particulier  de  mouvement  de  cible  radar.  Elle  est  applicable  a 
tout  mouvement  de  cible  realiste. 
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Executive  summary 


The  Inverse  Synthetic  Aperture  Radar  (ISAR)  is  a  technique  for  forming  high- 
resolution  images  by  utilizing  the  information  inherent  in  the  differential  target 
Doppler  that  results  from  rotation  of  an  object  (the  target)  relative  to  the  radar.  Its 
applications  include  imaging  and  discrimination  of  ships  and  aircraft  with  airborne 
or  land  based  radar  systems. 

Polarimetric  SAR/ISAR  imaging  is  an  effective  way  to  acquire  high  resolution  im¬ 
ages  of  targets  of  interest  at  long  range  and  as  such  is  an  irreplaceable  tool  in  the 
task  of  non-cooperative  target  recognition  of  both  ships  and  ground  moving  targets. 
One  of  the  significant  problems  with  SAR/ISAR  image  formation  of  a  moving  tar¬ 
get  is  the  assumption  of  time  invariance  of  the  Doppler  frequency  used  to  resolve 
the  image  in  the  cross  range.  Time  variant  Doppler  becomes  present  in  a  SAR  sig¬ 
nal  when  a  moving  target  is  maneuvering,  or  a  ship  is  pitching  and  rolling  during 
the  coherent  processing  interval  and  is  typically  referred  to  as  motion  error.  Us¬ 
ing  the  fast  Fourier  transform  to  cross  range  focus  an  image  with  this  motion  error 
present  will  cause  extensive  blurring  in  the  cross  range  and  leave  the  image  unre¬ 
cognizable  even  to  the  most  experienced  SAR  operators.  The  aim  of  this  report  is 
to  compensate  for  the  target  motion,  to  generate  an  image  equivalent  to  that  of  the 
stationary  target. 

The  Canadian  Air  Force  is  currently  upgrading  her  fleet  of  CP- 140  Aurora  mari¬ 
time  patrol  aircraft  to  possess  NCTR  through  ISAR  in  order  to  increase  the  capab¬ 
ility  of  the  Canadian  Forces  in  both  sovereignty  patrols  of  Canadian  territory  and 
the  protection  of  the  Canadian  Patrol  Frigates  and  allied  ships  operating  abroad  as 
a  coalition.  The  United  States  Navy’s  equivalent  aircraft,  the  P-3  Orion,  currently 
possesses  this  capability.  Therefore,  effective  ISAR  imaging  will  have  a  real  im¬ 
pact  in  the  decision  making  process  in  future  military  operations  involving  both 
Canadian  and  American  forces. 

In  this  report  the  adaptive  local  polynomial  Fourier  transform  is  employed  for  im¬ 
provement  of  the  ISAR  images  in  complex  reflector  geometry  cases,  as  well  as  in 
cases  of  fast  maneuvering  targets.  It  has  been  shown  that  this  simple  technique  can 
produce  significantly  improved  results  with  a  relatively  modest  calculation  burden. 
Two  forms  of  the  adaptive  LPFT  are  proposed.  The  adaptive  parameter  in  the  first 
form  is  calculated  for  each  radar  chirp.  The  additional  refinement  is  performed  by 
using  the  information  from  the  adjacent  chirps.  The  second  technique  is  based  on 
the  determination  of  the  adaptive  parameter  for  different  parts  of  the  radar  image. 
The  numerical  analysis  demonstrates  the  accuracy  of  the  proposed  techniques.  It  is 
important  to  note  that  the  proposed  technique  does  not  assume  any  particular  model 
of  radar  target  motion.  It  can  be  applied  for  any  realistic  motion  of  targets. 
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Sommaire 


La  technique  ISAR  (radar  a  synthese  d’ouverture  inverse)  permet  de  former  des 
images  a  haute  resolution  a  l’aide  de  1’ information  inherente  du  Doppler  differen- 
tiel  de  cible  que  produit  la  rotation  d’un  objet  (la  cible)  par  rapport  au  radar.  Ses 
applications  comprennent  1’imagerie  et  la  discrimination  de  navires  et  d’aeronefs 
au  moyen  de  systemes  radar  aeroportes  ou  bases  au  sol. 

Moyen  efficace  d’acquerir  des  images  a  haute  resolution  de  cibles  d’interet  a  longue 
distance,  l’imagerie  polarimetrique  SAR/ISAR  est  un  outil  irrempla?able  pour  la  re¬ 
connaissance  de  cibles  non  cooperatives  (NCTR),  y  compris  les  navires  et  les  cibles 
au  sol.  Un  probleme  important  de  la  formation  d’ images  de  cible  SAR/ISAR  tient 
a  ce  que  la  frequence  Doppler  utilisee  pour  resoudre  l’image  mobile  en  portee  late- 
rale  est  censee  etre  invariable  dans  le  temps.  La  frequence  Doppler  d’un  signal  SAR 
commence  a  varier  dans  le  temps  lorsqu’une  cible  mobile  effectue  des  manoeuvres, 
ou  en  presence  de  roulis  et  de  tangage  pendant  l’intervalle  de  traitement  coherent. 
On  parle  alors  generalement  d’"  erreur  de  mouvement  ".  L’utilisation  de  la  trans- 
formee  de  Fourier  rapide  pour  mettre  au  point  Fimage  en  portee  laterale  presentant 
cette  erreur  de  mouvement  produit  un  important  brouillage  de  la  portee  laterale  et 
rend  l’image  meconnaissable  meme  a  l’operateur  SAR  le  mieux  averti.  Le  present 
rapport  etudie  comment  compenser  le  mouvement  de  la  cible  afin  d’obtenir  une 
image  equivalant  a  celle  de  la  cible  immobile. 

La  Force  aerienne  du  Canada  est  en  train  de  mettre  a  niveau  la  flotte  d’avions  de 
patrouille  maritime  CP- 140  Aurora  en  les  dotant  de  la  fonction  NCTR  au  moyen  de 
1’ISAR  afin  d’augmenter  la  capacite  des  Forces  canadiennes  tant  dans  les  patrouilles 
de  protection  de  la  souverainete  du  territoire  canadien  que  dans  la  protection  des 
fregates  de  patrouille  canadiennes  et  des  navires  allies  faisant  partie  d’une  coalition 
a  l’etranger.  L’aeronef  equivalent  de  la  marine  americaine,  le  P-3  Orion,  est  deja 
dote  de  cette  capacite.  Par  consequent,  une  imagerie  ISAR  efficace  aura  une  forte 
incidence  sur  le  processus  de  prise  de  decisions  lors  d’operations  militaires  futures 
effectuees  par  des  forces  canadiennes  et  americaines. 

La  transformee  de  Fourier  polynomiale  locale  (TFPL)  adaptative  est  ici  utilisee 
pour  ameliorer  les  images  ISAR  dans  le  cas  de  reflecteurs  a  geometrie  complexe  et 
de  cibles  effectuant  des  manoeuvres  rapides.  II  a  ete  demonte  que  cette  technique 
simple  ameliore  sensiblement  les  resultats  au  prix  d’un  effort  de  calcul  relative- 
ment  modeste.  Deux  formes  de  la  TFPL  adaptative  sont  proposees.  Le  parametre 
adaptatif  de  la  premiere  forme  est  calcule  pour  chaque  impulsion  comprimee  radar. 
Le  resultat  est  raffine  a  l’aide  de  l’information  provenant  des  impulsions  compri- 
mees  adjacentes.  La  deuxieme  technique  se  fonde  sur  la  determination  du  parametre 
adaptatif  pour  differentes  parties  de  I’image  radar.  L’ analyse  numerique  montre 
1’exactitude  des  techniques  proposees.  Notons  que  la  technique  proposee  ne  sup- 
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1  Introduction 


The  inverse  synthetic  aperture  radar  (ISAR)  has  attracted  wide  interest  within  sci¬ 
entific  and  military  community.  Some  ISAR  applications  are  already  well  known 
and  studied.  However,  many  important  issues  remain  to  be  addressed.  For  example, 
suitable  enhancement  techniques  for  the  fast  maneuvering  radar  targets  or  targets 
with  fast  moving  parts  is  not  yet  known.  Also,  the  standard  approaches  based  on 
the  Fourier  transform  (FT)  fail  to  resolve  the  influence  of  closely  spaced  reflect¬ 
ors.  There  are  several  techniques  for  improvement  of  the  ISAR  radar  image  in  the 
case  of  fast  maneuvering  targets  or  in  the  case  of  objects  with  complex  reflector 
geometry.  Here  we  mention  only  two  groups  of  such  enhancement  techniques: 

•  techniques  that  adapt  transform  parameters  for  assumed  parametric  target  mo¬ 
tion  model  [1,6]; 

•  techniques  where  reflection  signal  components  are  parametrized,  while  the  sig¬ 
nal  components  caused  by  reflectors  are  estimated  by  using  some  of  the  well 
developed  parametric  spectral  estimation  tools  [2,  3]. 

Both  of  these  techniques  have  some  advantages,  but  they  also  have  some  drawbacks 
for  specific  applications.  The  first  group  of  techniques  is  strongly  based  on  the 
radar  target  geometry  with  an  assumed  motion  model.  These  techniques  could 
become  inaccurate  in  the  case  of  a  changing  motion  model.  The  second  group 
of  techniques  is  tested  on  simulated  examples.  However,  its  application  in  real 
scenarios,  where  signal  components  are  caused  by  numerous  scatterers,  could  be 
very  difficult.  Namely,  there  are  no  appropriate  methods  for  parameter  estimation 
of  signals  with  a  very  large  number  of  components. 

In  this  report  we  propose  a  modification  of  the  first  group  of  research  techniques. 
The  adaptive  local  polynomial  Fourier  transform  (LPFT)  is  used.  The  adaptive 
coefficients  are  calculated  for  each  considered  chirp  in  the  radar  signal  mixture. 
It  is  important  to  note  that  the  proposed  technique  does  not  assume  any  particu¬ 
lar  model  of  radar  target  motion.  The  adaptive  parameters  are  estimated  for  each 
scattering  point  independently.  Based  on  the  analysis  of  the  signal  obtained  from 
the  target  we  consider  some  simplifications  in  the  process  of  calculation  of  the  ad¬ 
aptive  transform.  In  this  way  we  keep  the  calculation  burden  to  within  reasonable 
limits.  The  two  techniques  for  the  enhancement  of  the  radar  image  by  using  the 
LPFT  are  considered.  The  first  one  is  based  on  the  information  obtained  from  each 
chirp  separately  and  on  the  possible  refinement  by  combining  results  from  various 
chirps.  The  second  technique  is  based  on  the  detection  of  regions-of-interest  in 
the  range/cross-range  plane  and  on  the  determination  of  the  optimal  LPFT  for  each 
detected  region  [19]. 
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The  report  is  organized  as  follows:  The  target  and  radar  signal  modeling  is  dis¬ 
cussed  in  Section  2.  The  proposed  methods  are  introduced  in  Section  3.  The  simu¬ 
lation  study  is  given  in  Section  4.  The  conclusions  are  given  in  Section  5. 
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2  Radar  Signal  Model 


Consider  a  radar  signal  consisting  of  M  continuous  wave  coherent  pulses: 


A/-1 

VA i(t)  =  v0 (t  -  mTr),  (1) 

m=0 

where  v0(t)  is  the  basic  impulse  limited  within  the  interval  —Tr/2  <  t  <  Trj 2. 
The  linear  frequency  modulated  (FM)  signal  is  used  in  our  simulations  as  a  basic 
impulse:  vo(t)  =  exp(jnBt2 /Tr),  where  B  is  the  bandwidth  control  parameter 
while  Tr  is  the  pulse  repetition  time.  The  alternative  radar  model  used  in  practice 
has  radar  pulses  with  stepped  frequencies,  however,  this  model  is  not  studied  in  this 
report.  The  defocusing  effect  considered  in  this  report  and  time- frequency  (TF) 
signatures  of  the  obtained  radar  signals  have  similar  behavior  for  these  two  forms 
of  radar  signals  [4,  5]. 

The  signal  emitted  toward  the  radar  target  can  be  written  as: 


u(t)  =  ej2^tvM(t),  .  (2) 

where  /0  is  the  radar  operating  frequency.  The  received  signal,  reflected  from  a 
single  reflector  target  at  distance  d(t),  is  delayed  for  2 d(t)/c,  with  a  being  the 
reflection  coefficient  and  c,  being  the  propagation  rate: 


uR(t)  =cru(t-  2d{t)/c) .  (3) 

The  demodulation  of  the  received  signal  can  be  performed  by  multiplying  the  re¬ 
ceived  with  the  transmitted  signal  u(t ): 

q(t)  =  au*(t  —  2  d(t)/c)u(t)  — 

A/-1  A/  — 1 

aexp(j4irf0d(t)/c)  v^t  -  2 d(t)/c-  mTr)  ^  ~  ™Tr  ~  T0).  (4) 

m.= 0  m= 0 

The  parameter  T0  is  used  in  radar  imaging  for  compensation  of  target  distance.  For 
properly  selected  T0  and  after  highpass  filtering,  the  signal  q(t)  can  be  approxim¬ 
ately  written  as: 
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A/-1 


M— 1 


q(t)  ~  aexp(j4irf0d{t)/c)  ^  v^t  -  2d(t) / c  -  mTr)v0(t  -  mTr)  =  ^  q(m,t), 


m= 0 


m— 0 


(5) 


where 


q(m,t )  =  crexp(j4:i:fQd(t)/c)vQ(t  —  2  d(t)/c  —  mTr)vo(t  —  mTr ), 
t  G  [(m  —  1/2 )Tr,  (rn  +  l/2)Tr), 

=  crexp(j47r/0d(f)/c)  exp(y'47r.Brf(f)(£  —  mTr)/(cTr ))  exp(—j7rB(2d(t)/c)2/Tr). 

(6) 

Keeping  in  mind  that  £?  C  /0.  we  can  neglect  exp(—jirB(2d(t)/c)2/Tr)  with  re¬ 
spect  to  the  other  two  components.  The  value  of  q(m,t )  can  approximately  be 
written  as: 


q{rn,  t)  «  a  exp(j4nf0d(t)/c )  exp(j4irBd(t)(t  -  mTr)/(cTr)).  (7) 
This  signal  is  commonly  given  in  the  form: 


q(m,  t)  «  a  exp(j4nf0d(T  +  m.Tr)/c. )  exp(/47rZ?fi(r  +  mTr)r/ (cTr)),  (8) 

where  t  =  r  +  mTr.  Parameter  r  G  [— Tr/2,  Tr/2)  is  referred  to  as  fast-time, 
while  m  =  0, 1, ....  M  —  1,  is  called  the  slow-time  coordinate.  Commonly,  in  actual 
radar  systems,  signals  are  discretized  in  fast-time  coordinate  with  sampling  rate 
Ts  —  Tr/N,  t  =  nTs,  where  n  G  [— JV/2,  N/2).  The  classical  radar  setup  assumes 
that  the  radar  target  position  is  a  linear  function  of  time  d(t)  =  D0  +  Vt.  Then  the 
radar  model  produces: 


q{m,  r)  «  cr  exp(j4irf0D0/c)  exp(j4irVmf0Tr/c)  exp(j4nTBD0/(cTr)).  (9) 
A  two-dimensional  (2D)  FT  of  this  signal  over  m  and  r  is  approximately: 

r  M—l 

Q(ujT,um)=  /  ?(m>  T)e~3“TT~JU}'nmdT  « 

m= 0 
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[  y  crexp(jinfc,D0/c)exp(j4i:VmJ0Tr/c)exp(jiTrTBD0/(cTr))e  ^ tT  mmdr 

^ T  m= 0 

„  A/-X 

p(j47r/0Z>o/c)  /  exp(j47TT SD0/ (cTr))e~JU)rTdr  y  exp(j4i:VmfQTr/c)e~3U” 

J  T  _ n 


=  a  ex 


/„  n  ,  n  /  \xt  a  on  rr^sin{(ujm-47rVf0Tr/c}M/2) 

=  {2ir)aexp{j47rf0D0/c)5(ujT  -  4~BD()/(eTr))  . . -  —  -  - 

sin  ( yUJfjj  4nvj()Tr/cj/2 ) 


xe 


—j(tAj,n—4nVfoTr/c)(M—l)/2 


(10) 


\Q(ujT,ojm)\  =  | (2tt)ct exp(j47r/o-D0/c)5(cur  -  47r.RDo/ (cTr)) 

x  S*I1((a;w  —  4nVf0Tr/c)M/2)  ?-(a>m-47rVfnTr/c)(A/-l)/2 
sin((wm  -  4-kV f0Tr/c)/2) 

sin((wm  -  4ttV f0Tr/c)M/2) 


sin((wm  -  AnV f0Tr/c)/ 2) 

\2ir\\a\\'exp(j4nfoD0/c)\\8(uT-4irBDo/(drr))\ 


sin((o;m  -  4nV f0Tr/c)/2) 


x  |e-j("m-4wV/0rr/c)(A/-l)/2| 

=  1(2 - 4,BPo/(cr,.))sin((^-4f  A7l/^^)| 

sm((wrn  -  AnV foTr/c)/2) 

«  (27r)<xA(wT  -  AvBDq/ {cT,))M8{ujm  -  4tt Vf0Tr/c) 


since 


sin(oM) 

sin(a) 


A/A(a) 


for  relatively  large  M. 

Therefore,  for  large  AT  we  can  write  the  magnitude  of  Q(uiT,uim)  as: 
|<5(wt,w„,)|  «  (2-)cr(5(wT.  -  4ttBD0/ (cTr ) ) M A (um  -  2 Vf0Tr/c). 


(11) 


(12) 


The  distance  can  be  approximately  written  as  d(t)  «  R(t)+xp  cos {9(t))+yp  sin (9(t)), 
where  R(t)  is  the  distance  of  the  target  rotation  center  from  the  radar,  where  co¬ 
ordinates  of  the  scatterer,  for  r  =  0,  are  {xp,  yp)  (see  Figure  1).  The  coordinate 
system  is  formed  in  such  a  way  that  the  coordinate  x  is  the  line  of  sight.  Assume 
constant  rotation  velocity  0(t)  =  uRt,  with  relatively  small  angular  movement  of 
the  target  \ujrTt\  <C  1  (it  implies  that  cos (0(f))  ~  1  and  sin(0(f))  ss  0).  Accord¬ 
ing  to  the  introduced  conditions  d(t)  «  xp  and  v(t)  =  d'(t )  =  — xp6'(t )  sin(0(£))  + 
yp0'(t )  cos {d(t))  ~  yp0'(t)  cos (9(t))  ~  Commonly,  it  is  assumed  that  R(t)  is 
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Figure  1:  Illustration  of  the  radar  target  geometry. 


compensated  by  adjusting  T0  in  (4).  Thus,  we  will  not  consider  it  in  our  algorithm. 
Then  \Q(ujT,ijjm)\  can  be  written  as: 


\Q{uT,um)\  «  (2iv)a MS (coT  -  47r Bxp/(cTr))5(um  -  4nypuiRf0Tr/c ) 

=  (2tt)<jM5(u)t  -  cixp)5(um  -  c2yp).  (13) 

It  represents  the  ISAR  image  of  scatterer  ( xp ,  yp)  for  a  given  instance  under  in¬ 
troduced  assumptions.  Note  that  the  constants  that  determine  the  resolution  of  the 
radar  image  are  given  by  Ci  =  4irB/(cTr)  and  c2  =  4iuxRf0Tr/c.  The  radar  image 
is  formed  as  a  superposition  of  radar  images  of  all  scatterers  (xp,  yp),  p  =  1,  2, ...,  P. 
It  is  approximately  given  as: 


p 

\Q(u)T,u]ni)\  =  y^(27r)cr p6(ujt  -  cixp)5(ujm  -  c2yp ),  (14) 

p=i 

where  op  is  the  reflection  coefficient  that  corresponds  to  the  p- th  scatterer  point. 

In  numerous  cases  we  cannot  assume  that  the  radar  model  can  be  simplified  in  the 
previously  described  manner.  For  example,  the  radar  target  can  be  very  fast,  or  the 
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model  of  radar  target  motion  can  be  more  complicated  (for  example  3D  motion). 
Then,  instead  of  complex  sinusoids  given  by  (9),  the  polynomial  phase  of  the  radar 
signal  model  can  be  written  as, 


q(m,T )  =  <Jp  exp  ,  (15) 

where  parameters  am,i  depend  on  the  considered  chirp  and  scatterer  motion.  For 
example,  for  the  target  motion  model  d(t)  =  D0  +  V0t  +  At2/ 2,  where  A  is  accel¬ 
eration  of  the  target,  the  coefficients  g,„5;  are  approximately  calculated  as  follows. 

Start  with  the  Equation  (8),  q(rn,  r)  ~  a  exp(j4Tr fod(r  +  mTr) / c)  exp(j4irBd(T  + 
mTr)T / (cTr)) ,  and  introduce  d(t)  =  D0  +  V0t  +  At2/ 2: 


q(m,T )  «  <7  exp (^^0  [D0  +  Vo (r  +  mTr)  +  A(t  +  mTr)2/ 2]) 

x  exp  (^^crp  ^  iDo  +  Ki(r  +  m-Tr )  +  A(t  +  mfr)2/2]rj  . 
Coefficients  amj  are  approximately  equal  to 

dm,  o  =  —[Do  +  VotnTr  +  d(mTr)2/2], 
c 

terms  with  r  are 

am,  i  =  ^[E0  +  AmTr]  +  i^-[D0  +  y0mTr  +  A(mTr)2/2], 

c  (cir) 


terms  with  r2  are 


terms  with  r3  are 


8tt  (  j,  A  V[ j  , 

dm,  2  =  —  (fo^  ■*"  ^ m  )  ’ 


dm,  3  — 


12itBA 
cTt  ’ 


(16) 


and  am., i  =  0  for  l  >  3.  Some  terms  of  these  coefficients  can  be  ignored,  but  in 
general  it  is  not  as  simple  as  in  the  case  when  we  can  assume  that  the  scatterer 
position  is  a  linear  function.  The  situation  becomes  even  more  difficult  in  the  case 
when  the  target  model  is  not  a  simple  rotating  model.  When  the  target  model  is 
complex,  a  very  complicated  relationship  between  the  position  of  scatterers  ( xp ,  yp) 
and  the  coefficients  of  the  polynomial  in  the  signal  phase  can  be  established.  Also, 
the  polynomial  that  should  be  used  to  accurately  estimate  the  signal  phase  will  be 


DRDC  Ottawa  TM  2006-185 


7 


of  very  high  order.  The  radar  image  obtained  by  using  the  2D  FT  of  the  signal 
with  higher  order  polynomial  becomes  spread  (defocused)  in  the  range/cross-range 
domain  (wr,  um).  The  goal  of  the  ISAR  signal  processing  is  to  obtain  the  focused 
radar  image,  i.e.,  to  remove  the  influence  of  the  higher  order  polynomial  in  signal 
phase  of  each  component. 

Usually,  it  is  assumed  that  the  modeling  of  coefficients  is  possible  based  on  the 
target  motion  model.  In  this  case,  instead  of  all  possible  parameters,  only  paramet¬ 
ers  of  the  motion  model  should  be  used  in  order  to  perform  an  enhancement  of  the 
radar  image. 

The  first  group  of  techniques  for  the  enhancement  of  radar  images  is  based  on  this 
concept.  One  such  approach  is  described  in  [6]  where  it  is  assumed  that  the  radar 
target  scatterers  can  be  modeled  with  a  relatively  simple  motion  model  which  as¬ 
sumes  that  the  velocity  increases  or  decreases  linearly  (or  that  angular  velocity 
changes  in  linear  manner)  within  the  repetition  time.  After  estimating  the  accel¬ 
eration  of  the  target,  the  variation  in  the  velocity  is  compensated  from  the  signal 
and  finally  a  focused  radar  image  is  obtained.  It  corresponds  to  removing  the  influ¬ 
ence  of  acceleration  from  (15).  However,  these  techniques  are  very  sensitive  to  any 
variations  from  the  assumed  motion  model.  They  cannot  be  used  for  3D  motion 
models. 

Alternative  techniques  are  based  on  estimation  of  all  the  amj  coefficients  in  the 
polynomial  phase  of  the  received  signal  [2,  3].  These  techniques  are  usually  based 
on  iteration,  removing  of  the  lower  order  coefficients  from  the  signal  phase  in  or¬ 
der  to  estimate  the  highest  order  coefficient.  Then,  the  estimation  of  lower  order 
coefficients  is  performed  by  using  the  same  procedure  but  for  the  dechirped  signal. 
It  means  that  the  error  in  the  estimation  of  the  highest  order  coefficient  propag¬ 
ates  toward  the  lower  order  coefficients.  Furthermore,  it  has  recently  been  shown 
that  these  procedures  are  biased  for  multicomponent  signals  and  that  the  dechirp¬ 
ing  procedure  used  to  produce  a  signal  suitable  for  estimation  of  the  lower  order 
coefficients  introduces  an  additional  source  of  errors  for  multicomponent  signals. 
These  techniques  are  also  time  consuming  and,  as  far  as  we  know,  never  applied 
to  signals  with  a  large  number  of  components.  Numerous  components  caused  by 
target  scatterers  could  appear  in  the  radar  signal. 

A  novel  technique  for  enhancing  radar  images,  that  introduces  just  one  new  adaptive 
parameter  in  the  FT  expression  for  each  received  signal,  is  introduced  in  the  next 
section.  For  each  chirp,  only  one  parameter  of  the  transform  should  be  estimated. 
The  second  important  property  of  this  technique  is  in  the  fact  that  we  do  not  assume 
any  particular  target  motion  model.  Therefore  it  can  be  applied  for  any  realistic 
motion  of  targets. 
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3  Adaptive  Local  Polynomial  FT 


In  this  section  we  introduce  the  LPFT  as  a  tool  for  the  autofocusing  of  ISAR  images. 
Two  forms  of  the  adaptive  LPFT  are  proposed.  The  first  form  can  be  applied  to 
each  chirp  component  separately  with  possible  refinements  by  using  information 
from  the  adjacent  chirps  (Section  3.1).  The  second  form  performs  evaluation  of  the 
adaptive  LPFT  for  each  detected  region-of-interest  in  the  radar  image  (Section  3.2). 

3.1  First  form:  Adaptive  LPFT  for  radar  signals 

In  order  to  develop  this  approach  we  will  go  through  several  typical  cases  of  signals, 
starting  from  a  very  simple  case  to  the  much  more  complicated  cases.  The  improve¬ 
ment  in  signal  components  concentration  (focusing  radar  image)  is  performed  by 
estimating  the  signal  parameters  without  assuming  any  particular  motion  model. 
This  is  a  quite  different  approach  compared  to  the  methods  with  predefined  motion 
model  or  to  the  methods  where  estimation  is  performed  on  each  parameter  arn  ;. 

3.1.1  Linear  FM  signal  case 

The  simplest  case  of  monocomponent  linear  FM  signal 


q(m,  r)  =  a  exp  (,j[am,o  +  a„hlT  +  a7„,2r2/2])  (17) 

is  considered  first.  In  this  case,  the  dependence  on  m  in  parameter  indices  will  be 
removed  for  the  sake  of  notation  brevity.  Then,  the  signal  can  be  written  as: 


q{m,  t)  =  o  exp  (j[a0  +  axr  +  a2r2 / 2]) 

For  the  analysis  of  these  kind  of  signals,  we  can  use  the  LPFT  [7],  [8]: 


(18) 


/OO 

q(m,  t)w(t)  exp(— jot2/2)  exp(—  jojTT)dr,  (19) 

-OO 

where  w(t)  is  a  window  function  of  the  width  Tw,  w(t)  =  0  for  |r|  >  Twj 2. 

The  LPFT  is  ideally  concentrated  along  the  instantaneous  frequency  for  a  =  a2: 


/OO 

w(t)  exp  (j[a0  -|-  o,\T  +  a2r2/2]) 

■OO 


X 


DRDC  Ottawa  TM  2006-185 


9 


(20) 


exp(— j£jTr  -  ja2T2/2)dr 

/CO 

w(t)  exp(—  j(uT  —  ai))d,T 

-CO 

=  aejaoW{ujT  -  ai) 

where  W(u>T)  =  FT{w(t)}.  Function  F(ojt,  m;  a2)  is  highly  concentrated  around 
U)T  =  a\ ,  since  the  FT  of  common  wide  window  functions  (rectangular,  Flamming, 
Hanning,  Gauss)  is  highly  concentrated  around  the  origin  (in  our  simulations  the 
window  width  is  equal  to  the  repetition  rate  Tw  =  Tr).  The  radar  image  can  be 
obtained  from  F(ojt,  m;  a2)  by  evaluating  the  ID  FT  along  the  m— coordinate  axis: 

A/-1 

Q(ojT,ojm\a2)  =  ^2  F(uT,m,;a2)e~:,ulmm.  (21) 

m= 0 


3.1.2  Higher  order  polynomial  FM  signal 

For  the  higher  order  polynomial  signal: 


q{m,,T)  =  a  exp  (]4>,„(t))  =  crexp(i</»(r))  (22) 

the  LPFT  can  be  written  as: 

/OO 

crexp(j0(r))u;(r)  exp(—  jar1  /2)  exp(—jojrr)dT 

■OO 

/OO 

exp(j0(O)  +  j4/{0)r  +  j^"(0)r2/2  +  j^"'(0)r3/3!  +  ... 

■OO 

+jd/n\Q)T7,/n\  +  ...  —  jar2 /2  —  jujTT)w{T)d,T  (23) 

For  0^(0)  =  0  for  n  >  2  we  obtain  a  highly  concentrated  LPFT  for  a  =  <f>"( 0): 

F(uit,  m;  4>"{0))  =  a  exp(jcf)(0))W(cuT  -  4>'{ 0))  (24) 

The  second  derivative  of  the  signal  phase  is  commonly  called  chirp-rate  parameter. 

In  the  case  when  higher  order  derivatives  are  non-zero,  the  LPFT  will  not  be  ideally 
concentrated  and  there  will  be  some  spread  in  the  frequency  domain  caused  by  the 
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FT  of  terms  exp(j0w(O)r3/3!  +  ...  +  j(j>^'\0)Tn/n\  +  ...).  The  LPFT  forms  that 
can  be  used  to  remove  effects  of  the  higher  order  derivatives  from  the  signal  phase 
are  introduced  in  [7],  [8].  These  techniques  are  computationally  demanding  and 
difficult  for  application  in  the  ISAR  imaging  in  the  real-time. 

An  alternative  technique,  the  order  adaptive  LPFT,  is  proposed  in  [9].  The  spectral 
width  of  the  signal’s  FT  is  used  as  an  indicator  of  the  polynomial  phase  order. 
Namely,  the  proper  order  and  parameters  of  the  LPFT  are  applied  if  its  width  in  the 
frequency  domain  is  close  to  the  width  of  considered  window  function  W ( u)T ). 

The  algorithm  for  the  order  adaptive  LPFT  determination  can  be  summarized  as 
follows: 

•  It  begins  with  the  ordinary  FT  calculation  (zero-order  LPFT)  in  the  first  step. 
If  the  width  of  this  transform  in  the  frequency  domain  is  equal  to  the  window 
width,  it  means  that  the  image  is  already  focused  and  there  is  no  need  for  the 
LPFT  order  increase.  Otherwise  go  to  the  next  step. 

•  Use  the  first  order  LPFT  form  considered  in  this  report,  eq.(19).  If  the  width  of 
this  transform  in  the  frequency  domain  is  equal  to  the  window  width,  it  means 
that  the  image  is  focused.  If  the  LPFT  still  have  some  spread  we  should  intro¬ 
duce  a  new  parameter,  j3,  into  the  transform  (next  coefficient  in  the  LPFT  phase 
will  be  —/3t3/ 3!)  and  repeat  the  operation. 

This  very  simple  idea  could  be  used  for  signals  with  one  or  at  most  few  components. 
In  complex  multicomponent  signal  cases,  more  sophisticated  techniques,  based  on 
the  concentration  measures,  will  be  introduced  in  the  next  section. 

3.1.3  Concentration  measure 

From  the  derivations  given  above,  it  can  be  concluded  that  for  a  known  chirp- 
rate  parameter  we  can  obtain  a  focused  radar  image  (highly  concentrated  time- 
frequency  (TF)  representation).  Also,  it  can  be  seen  that  the  ISAR  imaging  based 
on  the  LPFT  for  a  known  chirp-rate  parameter  is  slightly  more  demanding  than 
the  standard  ISAR  imaging  since,  in  addition  to  the  standard  procedure,  it  requires 
multiplication  with  the  term  exp(—  jar2  j2).  The  next  question  is  how  to  determine 
a  value  of  the  parameter  a  that  will  produce  highly  concentrated  images.  There  are 
several  methods  in  the  open  literature.  Here,  the  concentration  measures  will  be 
used  [10,  11,  12].  Before  we  propose  our  concentration  measure,  some  properties 
of  the  LPFT  will  be  reviewed.  The  LPFT  satisfies  the  energy  conservation  property: 


m;  q )F*(coTl  m;  a)dojT 
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r*oo  poo  /»oo 


-OO  t/  — oo  «/  — oo 


/oo 

q(m,,  Ta)w(ra)  exp(-jaT2a/2 )  exp(-jwTr0): 

-OO 

<f  (m,  Tb)w{rb)  exp(jar2/2)  exp (jujTTb)dTadrbdxjT 

/OO  /*oo 

/  q{m,  Ta)w(Ta)  exp(— jar„/2)  x 
-oo  J  —  OO 

f/*(m,  Tb)w(rb)  exp(jaT2b/2)S(Ta  -  Tb)dTadrb 

/OO 

|<jr(m,  r)|2u)2(r)dr. 

-OO 


(25) 


Consider  now  the  measure  \F(uT,  m;  a)|7dwT  for  7  — >  0.  Assume  that  F(wx,  m;  a) 

is  concentrated  in  a  narrow  region  around  the  origin  in  the  frequency  domain: 


|F(wt,  m:,a)\  =  0  forwT  >  0/2.  (26) 

Then,  we  obtain: 

/OO 

|F(u;r,m;  o)|7dwr  =  Q.  (27) 

-OO 


We  can  see  that  the  considered  measure  is  smaller  in  the  case  of  signals  concentrated 
in  narrower  intervals  in  the  TF  plane.  Therefore,  this  type  of  measure  can  be  used  to 
indicate  concentration  of  the  TF  representation.  In  a  realistic  scenario,  where  signal 
side  lobes  and  noise  exist  within  the  entire  interval,  this  measure  with  7  =  0  cannot 
be  used,  since  it  will  produce  an  approximately  constant  value.  In  order  to  handle 
this  issue,  we  can  use  0  <  7  <  2  instead  of  7  =  0.  As  a  good  empirical  value  in 
our  analysis  we  adopted  7  =  1.  Accurate  results  can  be  achieved  for  a  wider  region 
of  7  6  [0.5, 1.5], 

The  concentration  measure  based  on  the  above  analysis  can  be  written  as: 


H(m,  a;  7) 


_ 1 _ 

JZo\F(UJr,rn;a)\~’duT\ 


(28) 


Fligher  concentrated  signals  will  be  represented  by  a  higher  value  of  concentration 
measure  (28).  This  concentration  measure  has  been  proposed  [11],  where  it  is  ana¬ 
lyzed  in  detail  and  compared  with  other  concentration  measures.  This  concentration 
measure  produces  accurate  results  for  multicomponent  signals  as  well. 


12 


DRDC  Ottawa  TM  2006-185 


3.1.4  Estimation  of  the  chirp  rate  based  on  the  concentration 
measure 

Determination  of  the  optimal  chirp  rate  parameter,  a,  can  be  performed  by  a  direct 
search  of  the  assumed  set  of  a  values: 


Qopt(m )  =  arg max  H(m,  a;  7)  (29) 

a£A 


This  search  procedure  is  performed  over  the  parameter  space  A  =  [0,  amax],  where 
Qmax  is  the  chirp-rate  that  corresponds  to  the  TF  plane  diagonal.  cvmax  can  be  writ¬ 
ten  as  amax  =  2tt(1/2Ts)/(NTs/2)  =  2tt /(NT^),  where  1/2T,  is  the  maximal 
frequency  that  can  be  achieved  with  sampling  rate  Ts  within  repetition  time  Tr, 
Ts  =  Tr/N.  The  direct  search  over  a  single  parameter  is  currently  considered  as 
an  acceptable  computational  burden.  However,  in  the  case  when  calculation  time 
is  critical,  faster  procedures  should  be  used.  For  example,  in  the  case  of  mono¬ 
component  signals  embedded  in  moderate  noise,  the  LMS  style  algorithm  can  be 
employed.  The  optimal  value  of  the  chirp-rate  parameter  can  be  evaluated  as: 


ai+i{m)  =  a,(m)  -  p 


a,(m)  —  a,:_i(m) 


(30) 


where  [H(m,  a,:(m);  7)— #(m,  a;_1(m);7)]/[a7-(m)— c^-^m)]  is  used  to  estimate 
gradient  of  the  concentration  measure  and  f.i  is  the  predefined  step  size.  This  form  of 
the  algorithm  has  been  implemented  and  applied  for  TF  representations  in  [1 1],  A 
very  fast  (but  sensitive  to  noise  influence)  technique  for  estimation  of  the  chirp-rate 
parameters  has  been  proposed  in  [13]. 

3.1.5  Multicomponent  signals 

The  previously  described  procedures  for  determination  of  the  adaptive  chirp-rate 
parameter  can  be  applied  when  the  reflected  chirp  can  be  represented  as  a  mono¬ 
component  FM  signal.  Furthermore,  the  same  procedure  can  be  applied  for  mul¬ 
ticomponent  signals  with  the  same  or  similar  second  derivatives  of  the  signal  phase, 
since  the  search  for  just  one  chirp-rate  parameter  needs  to  be  performed.  This  situ¬ 
ation  corresponds  to  closely  spaced  scatterer  points  in  the  radar  image  with  similar 
motion  trajectories. 

However,  a  modification  is  required  in  the  case  of  several  components,  with  dif¬ 
ferent  chirp-rates.  In  this  case  the  previously  described  algorithm  would  produce  a 
highly  concentrated  dominant  signal  component,  while  the  remaining  components 
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would  be  spread  across  the  TF  plane.  The  method  proposed  in  [14]  is  based  on  the 
calculation  of  an  adaptive  transform,  as  a  weighted  sum  of  the  LPFTs: 

1 

m,  a; 7 )da 


*  AD\wt,  ni) 


rx 


/: 


F(cjt,  m;  a)H(m,  a;  ^)da  (31) 


where  the  weighted  coefficients  are  proportional  to  the  concentration  measure.  This 
method  had  produced  good  results  for  signals  with  components  of  similar  mag¬ 
nitudes.  Flowever,  if  signal  components  significantly  differ  in  amplitude,  the  results 
are  not  satisfactory.  Namely,  signal  components  with  smaller  amplitude  would  be 
additionally  attenuated.  In  order  to  avoid  this  drawback,  we  will  use  the  following 
adaptive  local  polynomial  FT: 


p 

Fad{lot,  m)  =  F(uit,  m;  af (m))  (32) 

i= 1 

where  the  first  adaptive  frequency  is  estimated  as: 


c*i  (m)  =  arg  max  ( rn ,  a;  7) 

a 


(33) 


with  H^(rn,  a;  7)  =  H(m,a;  7),  given  with  (28)  and  set  i  —  1.  After  detection 
of  the  first  component’s  chirp-rate,  values  of  H(m,  a ;  7)  in  a  narrow  zone  around 
07(772 )  are  neglected,  and  the  search  for  the  next  maximum  is  performed.  Each 
iteration  in  this  procedure  could  be  described  into  two  steps: 


H^\m,  o;  7)  '= 


^(m,o;7)  |o  —  a7;(?n)|  >  A 


0  otherwise 

a,+i(m )  =  arg  max  H^im,  o;  7),  i  =  i  +  1. 


(34) 


This  procedure  should  be  stopped  after  the  maximal  value  of  arg  maxa  (m,  o;  7) 
becomes  smaller  than  an  assumed  threshold.  We  set  that  the  threshold  is  25%  of 
max„  H^(m,  a;  7),  i.e.,  25%  of  concentration  measure  before  we  start  with  peel¬ 
ing  off  components.  Note  that  the  parameter  A  should  be  selected  carefully  so 
that  the  next  recognized  component  is  not  just  a  “side  lobe”  of  the  previous  strong 
component.  In  the  case  when  components  have  chirp-rates  close  to  each  other,  it  is 
enough  to  recognize  a  single  chirp-rate,  since  the  proposed  approach  will  improve 
concentration  of  all  the  components  with  similar  chirp-rates.  In  our  experiments  we 
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assumed  that  the  number  of  components  with  different  chirp-rates  for  considered 
radar  chirp  cannot  be  larger  than  8  and  we  selected  that  A  =  amax/16  =  7r/(8NTA2) 
(see  also  section  3.1.4).  It  produces  accurate  results  in  all  of  our  experiments.  Note 
that  an  alternative  method  for  evaluation  of  the  LPFT  is  proposed  in  [15]. 

3.1 .6  Combination  of  the  results  from  various  radar  chirps 

In  the  case  of  radar  signals  we  can  assume  that  scatterers  closely  spaced  in  the 
range/cross-range  plane  have  similar  motion  parameters.  This  means  that  for  chirps 
with  similar  chirp  number  we  can  take  similar  values  of  the  chirp-rate  parameter. 
The  chirp-rate  estimated  for  the  m-th  chirp  can  be  used  with  a  small  error  for  the 
next  chirp  signal,  without  recalculating  the  concentration  measure.  This  simpli¬ 
fied  technique  was  accurate  for  simple  simulated  reflector  geometry.  In  the  case 
of  complex  reflector  geometry,  with  numerous  closely  spaced  components,  inac¬ 
curate  chirp-rate  parameter  estimates  are  obtained  for  several  chirps.  Usage  of  one 
chirp-rate  for  the  following  chirps  causes  a  propagation  effect  error.  Therefore,  the 
concentration  measure  must  be  calculated  and  chirp-rate  parameter  should  be  es¬ 
timated  for  each  chirp.  In  order  to  refine  the  results  further,  non-linear  filtering  of 
the  obtained  chirp-rates  is  performed.  Assuming  that  the  chirp-rate  parameter  a(m) 
is  estimated  for  each  chirp,  the  nonlinear  median  filter  can  be  calculated  as: 


a(rn)  —  median{o(m  4-  i),  i  G  [— r,  r]}  (35) 

where  2 r  +  1  is  the  width  of  the  used  median  filter.  Note  that  other  filters  with  the 
ability  to  remove  impulse  noise  can  be  used  here  instead  of  the  median  filter,  for 
example  the  a-trimmed  mean  filters  [16,  17]. 


3.2  Second  form:  Adaptive  LPFT  for  regions  of  the 
radar  image 

Methods  for  adaptive  calculation  of  the  radar  image  described  so  far  propose  evalu¬ 
ation  of  the  adaptive  parameter  for  each  considered  chirp  and  a  possible  refinement 
by  combining  results  obtained  on  closely  spaced  scatterers.  The  implicit  assump¬ 
tion  was  that  the  close  points  in  the  range/cross-range  domain  have  similar  chirp- 
rate  parameters.  In  order  to  have  a  more  robust  technique  that  is  able  to  deal  with 
more  challenging  motion  models,  we  propose  an  alternative  form  of  the  adaptive 
LPFT  with  2D  optimization  of  chirp  parameters.  In  defining  this  procedure,  we 
keep  in  mind  that  relatively  small  portion  of  the  radar  image  is  related  to  the  target. 
We  will  consider  just  the  part  of  the  radar  image  that  exceeds  a  threshold: 
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Ie  (c 0T ,  U)nl ) 


1  \Q(uT,ujm)\  >  emax{|Q(wr,wm)|} 


0 


otherwise. 


(36) 


The  region  L( uT,  um)  can  be  separated  into  non-overlapping  regions: 


Pz 

4(ut:  ^m)  ^m)  (37) 

i= 1 

where  I,(ujT.  u]m  )  n  I3  (ojt ,  ojm )  =  0  for  i  ^  j.  Note  that  the  number  of  separated 
regions  pE  depends  on  selected  threshold  e.  By  using  the  inverse  2D  FT  we  can 
calculate  signals  associated  with  the  region  /,( uT,u>m) 


qi(m,T)  =  IFT{Q(ujT,u)m)Ii(u>T,ujTr,)},  i  =  1,2,  ...,pe.  (38) 

Now,  we  can  assume  that  signal  qi(m,  r)  is  generated  by  a  single  reflector.  Then, 
we  can  perform  optimization  of  each  signal  <7,:(m,r).  Since  this  signal  is  already 
localized  in  the  range/cross-range  domain,  we  will  not  perform  optimization  for 
each  r  or  m,  but  only  optimization  with  a  single  chirp  function  for  each  region 

7»(wr,  ^ m ) • 


poo  M _ ) 

Fi{ur,uim;  oti )  =  T)  exp(-jQ,;r2/2  -  juTr  -  jumm)dr  (39) 

m= 0 


where 


ctj  =  arg  max  ■ 


°  IZo  Em=0  I Fi(uT,  Urn,  a)\~'duJT 


(40) 


The  radar  image  is  calculated  as  a  sum  of  the  adaptive  LPFT  F,:(wr,  ujm;  a,): 


Pz 

Fs^ADif^Ti^in)  ^  ^  F'/ (l UT .OJm,  OLj ) ■  (41) 

i=l 

In  our  experiments  we  obtain  very  good  results  for  e  over  a  relatively  wide  range 
for  numerous  radar  images. 
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However,  additional  optimization  can  be  done  based  on  the  threshold  e.  Here,  a 
three-step  technique  for  threshold  selection  is  considered.  In  the  first  stage  we 
consider  various  thresholds  e  e  E  and  calculate  F^ad{^ti  wm)  for  each  threshold 
from  the  set.  Then,  we  calculate  the  optimal  LPFT  as  F£>AD(uT,wm)  that  achieves 
the  best  concentration  over  e  €  E.  Since,  by  introducing  the  threshold  value,  we 
remove  a  part  of  the  range/cross-range  plane  (see  (36))  the  energy  of  F^ad^t,  win) 
should  be  normalized  to  the  energy  of  signal  above  the  specific  threshold: 


Fe^AD^^ti  m ) 


/“ o  E!_0  | Q{uT,Um)\2I£{Uu,Um)dut 


e  —  arg  max  — — - n— ; - 

6  f°°  V'AZ-I  I  rv 


-'s  ns 


In  this  procedure  the  transforms  F£i/i£)(a;T,a;TO),  e  G  H,  are  compared  under  un¬ 
equal  conditions  since  they  are  obtained  with  various  thresholds  e  and  they  could 
have  different  number  of  recognized  components.  The  obtained  adaptive  transform 
Fs<ad(ut,  uw)  could  be  poorly  concentrated  than  a  particular  FL  AD{uT,um)  from 
the  considered  set  of  e  values.  However,  this  radar  image  is  near-optimal  and  a 
small  additional  manual  adaptation  around  the  estimated  e  could  be  performed  in 
the  third  stage  of  this  procedure  to  obtain  the  optimal  image.  In  our  experiments 
the  obtained  i  is  underestimated  (see  Example  7).  Thus,  an  additional  search  could 
be  performed  over  higher  values  of  e. 
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4  Results 


Several  numerical  examples  will  be  presented  here  to  justify  the  presented  ap¬ 
proach.  Examples  1-4  are  generic  signals  representing  one  received  radar  chirp 
that  prove  that  the  adaptive  LPFT  can  be  used  to  produce  highly  concentrated  TF 
representation  for  the  following  ID  signals:  linear  FM,  sinusoidal  FM,  multicom¬ 
ponent  signal  with  similar  chirp-rates  and  multicomponent  signals  with  differing 
chirp-rates.  Examples  5  and  6  demonstrate  that  the  adaptive  LPFT  optimized  for 
each  chirp  signal  with  filtering  data  produced  by  adjacent  radar  chirps  gives  accur¬ 
ate  results.  Example  7  illustrates  the  second  adaptive  LPFT  algorithm  with  optim¬ 
ization  over  regions  of  the  radar  image. 


4.1  Example  1 

The  first  signal  that  will  be  considered  is  a  linear  FM  signal  f(t )  =  exp(j647rf2/2) 
embedded  in  Gaussian  noise  with  variance  cr2  =  1.  The  signal  is  sampled  with 
At  =  l/128sec.  The  FT  of  the  windowed  signal  with  a  Hanning  window  of  the 
width  T  =  2  sec  is  shown  in  Figure  2a.  It  can  be  seen  that  the  FT  is  spread.  Thus,  if 
this  signal  is  a  part  of  the  received  signals  reflected  from  a  target,  we  will  obtain  a 
defocused  radar  image.  Results  obtained  with  narrower  Hanning  windows  are  given 
in  Figure  2b.  The  improvement  could  be  observed  from  this  figure,  but  generally 
speaking  it  is  slight.  The  concentration  measure  (28)  for  7  =  1  is  presented  in 
Figure  2c,  with  marked  detected  chirp-rate  parameter.  Finally,  the  adaptive  LPFT  is 
given  in  Figure  2d  calculated  for  parameter  a  for  which  the  concentration  measure 
given  in  Figure  2c  is  maximized.  The  significant  improvement  achieved  by  the 
LPFT  is  obvious. 

4.2  Example  2 

The  second  signal  is  a  more  complex  sinusoidal  FM  signal:  f(t )  =  exp(jl6  sin(27rt)) 
The  signal  sampling  and  noise  environment  are  the  same  as  in  Example  1.  The  FTs 
with  a  wide  and  a  narrow  window  around  a  given  time  instant  (STFT),  [18],  are  de¬ 
picted  in  Figures  3a, b.  This  STFT  illustration  for  a  fixed  instant  corresponds  to  the 
radar  image  for  considered  m.  It  can  be  used  to  estimate  the  radar  image  depending 
on  different  chirp-rates.  Again  we  can  see  that  for  each  instance  this  representa¬ 
tion  is  spread  in  frequency  domain.  It  means  that  the  radar  image  obtained  based 
on  the  FT  for  signals  of  this  form  will  be  defocused.  The  adaptive  LPFT  with  a 
single  chirp-rate,  calculated  for  each  instant,  is  given  in  Figure  3c.  A  significant 
improvement  is  achieved.  Also,  it  can  be  noticed  that  the  representation  is  not  ideal 
in  the  region  with  higher  order  derivatives.  These  derivatives  can  be  removed  by 
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Figure  2:  Spectral  analysis  of  the  linear  FM  signal:  (a)  FT  with  a  wide  window;  (b) 
FT  with  a  narrow  window;  (c)  Concentration  measure;  (d)  Adaptive  LPFT. 
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Figure  3:  Time-frequency  analysis  of  the  sinusoidal  FM  signal:  (a)  STFT  with  a 
wide  window;  (b)  STFT  with  a  narrow  window;  (c)  Adaptive  LPFT;  (d)  Adaptive 
chirp-rate  parameter. 
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employing  higher  order  LPFT  form  [7] -[9].  The  adaptive  chirp  rate  is  given  in 
Figure  3d. 


4.3  Example  3 


A  three-component  signal:  /(f)  =  exp(/227rf2+j487rf)-l-exp(/327rf2)+exp(/427rf2— 
jA8nt)  is  considered  next.  The  STFT  with  a  wide  and  a  narrow  window  is  given 
in  Figure  4a, b.  The  adaptive  LPFT  calculated  as  in  the  case  of  monocomponent 
signal  is  given  in  Figure  4c.  It  can  be  seen  that  the  concentration  is  improved  for 
all  three  components.  The  component  in  the  middle  is  enhanced  the  best,  but  other 
components  with  similar  chirp  rates  are  also  improved.  The  adaptive  parameter  is 
given  in  Figure  4d.  This  case  corresponds  to  a  signal  obtained  from  several  scatter- 
ers  in  the  same  cross-range  with  similar  chirp-rates.  The  difference  in  chirp-rates 
of  these  components  in  fact  is  not  so  small,  it  is  30%  of  the  chirp-rate  of  the  middle 
component.  In  real-time  applications  it  represents  a  realistic  scenario  of  the  scat¬ 
tered  in  the  radar  image.  We  can  see  that  the  concentration  of  all  components  is 
satisfactory.  It  can  also  be  seen  that  accuracy  of  this  procedure  is  not  affected  by 
the  distance  between  scatterers.  The  same  accuracy  is  achieved  for  the  left  part  of 
Figure  4c,  where  we  assume  that  scatterers  are  far  from  each  other,  as  well  as  in 
the  right  part  of  this  illustration,  where  it  can  be  assumed  that  scatterers  are  close  to 
each  other. 


4.4  Example  4 


A  three-component  signal:  /(f)  =  exp(/ll7rf2+/487rf)+exp(/327rf2)-|-exp(j677rf2— 
j'487rf)  is  considered.  However,  in  this  case  the  chirp-rates  of  components  are 
quite  different  (difference  between  chirp-rates  is  more  than  60%  of  chirp-rate  of 
the  middle  component).  The  STFT  is  given  in  Figure  5a,  while  the  “adaptive” 
transform,  assuming  the  signal  has  a  single  chirp-rate,  is  given  in  Figure  5b.  It  can 
be  seen  that  in  each  instant,  the  transform  is  adjusted  to  one  component,  while  other 
components  remain  spread.  For  f  <  0.3  the  LPFT  is  highly  concentrated  for  the 
middle  component,  but  when  components  are  close  to  each  other  (corresponding 
to  closely  spaced  scatterers)  the  adaptive  chirp-rate  switches  several  times  between 
components.  The  adaptive  weighted  LPFT  (32)  is  given  in  Figure  5c.  It  can  be 
seen  that  all  components  have  improved  concentration  and  that  concentration  is  not 
influenced  by  the  distance  between  scatterers.  The  detected  adaptive  chirp-rates  are 
given  in  Figure  5d. 


DRDC  Ottawa  TM  2006-185 


21 


300 


Figure  4:  Time-frequency  analysis  of  the  multicomponent  signal:  (a)  STFT  with  a 
wide  window;  (b)  STFT  with  a  narrow  window;  (c)  Adaptive  LPFT;  (d)  Adaptive 
chirp-rate  parameter. 
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Figure  6:  Simulated  radar  image:  (a)  Results  obtained  by  the  FT;  (b)  Adaptive 
chirp-rate  parameter  as  function  of  m  (thick  line  is  linear  approximation);  (c)  Radar 
image  based  on  the  adaptive  LPFT. 
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4.5  Example  5 


The  simulated  radar  target  setup  according  to  the  experiment  in  [4]  is  considered. 
The  reflectors  are  at  the  positions  (x,y)  =  {(—2.5,1.44),  (0,1.44),  (2.5,1.44), 
(1.25,  -0.72),  (0,2.88),  (—1.25,0.72)}  in  meters.  The  high  resolution  radar  op¬ 
erates  at  the  frequency  /0  =  10.1GHz,  with  a  bandwidth  of  linear  FM  chirps 
B  =  300MHz  and  pulse  chirp  repetition  time  Tr  =  15.6ms.  The  target  is  at 
2km  distance  from  the  radar,  and  rotates  at  uR  =  4°/  sec .  The  nonlinear  rota¬ 
tion  with  frequency  Q  =  0.5Hz  and  amplitude  A  =  1.25°/  sec  is  superimposed, 
uiR(t)  =  ojr  +  v4sin(27rflt).  The  FT  based  image  of  the  radar  target  is  depicted 
in  Figure  6a.  The  radar  image  obtained  by  using  the  adaptive  LPFT  calculated  for 
each  chirp  separately  is  presented  in  Figure  6c,  while  the  adaptive  parameter  for 
each  chirp-signal  is  given  in  Figure  6b.  It  can  be  seen  that  the  adaptive  parameter 
linearly  varies  between  the  limits  of  the  target.  Flowever,  the  impulse  like  errors  in 
the  estimation  of  the  chirp-rate  can  be  observed  from  Figure  6b.  It  suggests  that  the 
improvement  of  the  results  can  be  achieved  by  filtering  chirp-rate  parameters. 


4.6  Example  6 

In  this  example  we  consider  simulated  Boeing-727  radar  data.  The  FT  based  image 
is  presented  in  Figure  7a.  It  can  be  seen  that  the  radar  image  is  defocused.  How¬ 
ever,  the  radar  imaging  based  on  the  adaptive  LPFT  determined  for  each  radar  chirp 
produces  a  significant  improvement  in  the  signal  representation,  Figure  7b.  In  order 
to  obtain  better  results  for  closely  spaced  reflectors,  we  consider  the  adaptive  chirp- 
rate  parameter  depicted  in  Figure  7c  as  a  dotted  line.  We  expected  that  removing 
impulse  like  disturbances  will  produce  better  results.  To  this  aim,  the  median  filter¬ 
ing  of  the  adaptive  parameter  is  performed.  In  order  to  evaluate  the  outcome,  the 
linear  interpolation  of  estimated  chirp-rates  is  performed  (The  linear  interpolation 
is  depicted  with  thick  line  in  Figure  7c).  The  result  obtained  with  these  parameters 
is  depicted  in  Figure  7d.  It  is  better  than  its  counterpart  in  Figure  7b  except  for 
the  nose  reflectors.  A  possible  reason  is  the  fact  that  the  received  signal  corres¬ 
ponding  to  these  scatterers  can  have  higher  order  polynomials  in  the  signal  phase. 
The  higher  order  LPFT  forms  [7]-[9]  could  be  used  for  these  scatterers  (see  Section 
3.1.2). 

4.7  Example  7 

In  this  example  we  consider  the  same  target  as  in  Example  5.  The  main  difference 
in  this  example  is  the  complex  motion  pattern  that  cannot  be  modeled  with  just  a 
rotation.  The  radar  image  calculated  by  using  the  2D  FT  is  presented  in  Figure 
8a.  Region-of-interest  Ie{ujT,ujm )  is  determined  by  (36)  with  the  threshold  set  to 
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Figure  7:  B727  radar  image:  (a)  Results  obtained  by  the  FT  based  method;  (b) 
Adaptive  LPFT  method;  (c)  Adaptive  chirp-rate  -  dotted  line;  Filtered  adaptive 
chirp-rate  -  dashed  line;  Linear  interpolation  of  filtered  data  -  solid  line;  (d)  Adaptive 
LPFT  with  interpolated  data. 
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Figure  8:  Simulated  radar  image  with  complicated  motion  pattern:  (a)  Results 
obtained  by  the  FT;  (b)  Regions  of  interest  Ie=o.os(ut,  wm)  with  three  recognized 
separated  regions;  (c)  Adaptive  LPFT  based  on  region  optimization  with  e  =  0.05, 
^e=o.o5 (wt )  um);  (d)  Regions  of  interest  4=o.2o(wt,  com)  with  six  recognized 
separated  regions;  (e)  Adaptive  LPFT  based  on  region  optimization  with  e  =  0.20, 

F r=0.20 y^t i  kJ»n)- 
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Figure  9:  Adaptive  LPFT  with  adaptive  threshold:  (a)  Concentration  measure  for 
various  threshold  levels.  Optimal  threshold  value  is  depicted  with  dashed  line,  (b) 
Adaptive  LPFT  with  adaptive  threshold. 


e  =  0.05.  Three  separated  regions  are  detected  in  the  radar  image,  denoted  in 
Figure  8b,  in  different  shades  of  gray.  The  region  denoted  with  1  corresponds  to 
three  radar  scatterers.  Since  these  three  scatterers  move  in  a  similar  manner,  the 
concentration  of  these  components  is  significantly  improved  (see  Figure  8c)  with 
respect  to  the  radar  image  calculated  with  2D  FT.  The  region  denoted  with  num¬ 
ber  3  corresponds  to  two  radar  scatterers.  In  this  case,  the  concentration  of  one  of 
the  components  from  the  region  is  improved,  while  the  other  component  remains 
spread.  The  reason  is  in  fact  that  these  scatterers  move  in  a  quite  different  manner. 
When  we  apply  threshold  e  =  0.2,  we  obtain  6  regions  of  interest  that  correspond 
to  6  radar  scatterers  (Figure  8d).  The  resulting  radar  image  is  focused  for  all  scat¬ 
terers  (Figure  8e).  The  threshold  e  could  be  set  in  an  empirical  manner.  However, 
a  procedure  for  threshold  optimization  could  be  very  helpful.  The  concentration 
measure  of  the  adaptive  LPFT  for  various  threshold  levels  is  depicted  in  Figure  9a 
and  the  obtained  value  in  the  optimization  procedure  is  e  «  0.155.  The  LPFT  form 
with  adaptive  threshold  is  shown  in  Figure  9b.  It  can  be  seen  that  the  radar  im¬ 
age  obtained  in  Figure  9b  is  slightly  worse  than  the  radar  image  with  additionally 
adjusted  threshold  Figure  8e. 
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5  Conclusion 


The  adaptive  local  polynomial  Fourier  transform  based  method  for  enhancement  of 
defocused  radar  images  has  been  proposed.  The  adaptive  parameters  in  the  trans¬ 
form  are  obtained  by  using  a  simple  concentration  measure.  For  monocomponent 
and  multicomponent  signals  with  similar  chirp-rates,  a  single  chirp-rate  parameter 
is  estimated  for  each  chirp.  For  multicomponent  signals  with  different  chirp-rates, 
an  adaptive  weighted  local  polynomial  FT  should  be  employed.  It  has  been  shown 
that  the  ISAR  images  could  be  improved  by  combining  results  achieved  from  vari¬ 
ous  chirps.  For  targets  with  very  complex  motion  pattern,  the  separation  of  the  radar 
image  in  regions-of-interests  and  optimization  of  the  radar  signal  within  regions  is 
proposed.  The  proposed  technique  does  not  assume  any  particular  model  of  radar 
target  motion.  It  can  be  applied  for  any  realistic  motion  of  targets. 


DRDC  Ottawa  TM  2006-185 


29 


References 


1.  Y.  Wang,  H.  Ling,  V.  C.  Chen:  “ISAR  motion  compensation  via  adaptive  joint 
time-frequency  techniques,”  IEEE  Trans.  Aer.  Ele.  Sys.,  Vol.  34,  No.  2,  Apr. 
1998,  pp.  670-677. 

2.  S.  Barbarossa,  A.  Scaglione,  G.  B.  Giannakis:  “Product  high-order  ambiguity 
function  for  multicomponent  polynomial-phase  signal  modeling,”  IEEE  Trans. 
Sig.  Proc.,  Vol.  46,  No.  3,  Mar.  1998,  pp.  691-708. 

3.  A.  Quinquis,  C.  Ioana,  E.  Radoi:  “Polynomial  phase  signal  modeling  using 
warping-based  order  reduction,”  in  Proc.  of  ICASSP ' 04 ,  Vol.  2,  May  2004,  pp. 
741-744. 

4.  S.  Wong,  E.  Riseborough,  and  G.  Duff:  “Experimental  investigations  on  the 
distortion  of  ISAR  images  using  different  radar  waveforms,”  Defence  R&D 
Canada  -  Ottawa,  DRDC  Ottawa  TM  2003-196,  2003. 

5.  S.  Wong,  E.  Riseborough,  and  G.  Duff:  “Distortion  in  the  ISAR  (inverse 
synthetic  aperture  radar)  images  from  moving  targets,”  in  Proc.  of  IEEE 
ICIP' 2004,  Vol.  I,  pp.  25-28,  2004. 

6.  T.  Thayaparan,  G.  Lampropouols,  S.  K.  Wong  and  E.  Riseborough, 
“Application  of  adaptive  joint  time-frequency  algorithm  for  focusing  distorted 
ISAR  images  from  simulated  and  measured  radar  data,”  IEE  Proc.  Radar 
Sonar  Navig.,  Vol.  150,  No.  4,  Aug.  2003,  pp.  213-220. 

7.  V.  Katkovnik,  “A  new  form  of  the  Fourier  transform  for  time-frequency 
estimation,”  Sig.  Proc.,  Vol.  47,  No.  2,  pp.  187-200,  1995. 

8.  V.  Katkovnik,  “Local  polynomial  periodogram  for  time-varying  frequency 
estimation,”  South  Afr.  Stat.  Jour.,  Vol.  29,  No.  2,  pp.  16,  168-195. 

9.  LJ.  Stankovic,  S.  Djukanovic,  “Order  adaptive  local  polynimial  FT  based 
interference  rejection  in  spread  spectrum  communicaton  systems,”  in  Proc.  of 
IEEE  WISP  2003. 

10.  R.  G.  Baraniuk,  P.  Flandrin,  A.  J.  E.  M.  Jensen,  O.  J.  J.  Michel,  “Measuring 
time-frequency  information  content  using  Renyi  entropy,”  IEEE  Trans.  Inf. 

Th.,  vol.  47,  no.  4,  May  2001,  pp.  1391-1409. 

11.  LJ.  Stankovic,  “A  measure  of  some  time-frequency  distributions 
concentration,”  Sig.  Proc.,  vol.  81,  no.  3,  Mar.  2001,  pp.  621-631. 

12.  T.  H.  Sang,  W.  J.  Williams,  “Renyi  entropy  and  signal  dependent  optimal 
kernel  design,”  in  Proc.  ICASSP,  vol.  2,  1995,  pp.  997-1000. 


30 


DRDC  Ottawa  TM  2006-185 


13. 1.  Djurovic,  LJ.  Stankovic:  “Moments  of  multidimensional  polynomial  FT,” 
IEEE  Sig.  Proc.  Let. ,  Vol.  1 1,  No.  1 1,  Nov.  2004,  pp.879-882. 

14.  M.  Dakovic,  I.  Djurovic,  LJ.  Stankovic,  “Adaptive  local  Fourier  transform”,  in 
Proc.  of  EUSIPCO’2002,  Toulouse,  France,  Vol. II,  pp. 603-606. 

15.  Y.  Wei,  G.  Bi,  “Efficient  analysis  of  time-varying  multi-component  signals 
with  LPTFT,”  Jour.  Appl.  Sig.  Proc.,  No.  8,  2005,  pp.  1261-1268. 

16. 1.  Pitas,  A.  N.  Venetsanopoulos,  Nonlinear  digital  filters:  Principles  and 
applications,  Kluwer  Academic,  1990. 

17. 1.  Djurovic,  LJ.  Stankovic,  J.  F.  Bohme,  “Robust  L-estimation  based  forms  of 
signal  transforms  and  time-frequency  representations,”  IEEE  Trans.  Signal 
Processing,  Vol.  51,  No.  7,  July  2003,  pp.  1753-1761. 

18.  J.  B.  Allen  and  L.  R.  Rabiner,  “A  unified  approach  to  short-time  Fourier 

analysis  and  synthesis,”  Proc.  IEEE,  vol.  65,  no.  11,  pp.  1558-1564,  Nov.  1977. 

19. 1.  Djurovic,  T.  Thayaparan,  and  LJ.  Stankovic,  "Adaptive  local  polynomial 
Fourier  transform  in  ISAR,  J.  on  Applied  Signal  Processing,  in  press,  2006. 


DRDC  Ottawa  TM  2006-185 


31 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 
(highest  classification  of  Title,  Abstract,  Keywords) 

DOCUMENT  CONTROL  DATA 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  document  is  classified) 


3.  TITLE  (the  complete  document  title  as  indicated  on  the  title  page.  Its  classification  should  be  indicated  by  the  appropriate 
abbreviation  (S,C  or  U)  in  parentheses  after  the  title.) 

Focusing  ISAR  Images  using  the  Adaptive  Local  Polynomial  Fourier 
Transform  (U) 

4.  AUTHORS  (Last  name,  first  name,  middle  initial) 

Thayaparan,  Thayananthan 


6a.  NO.  OF  PAGES  (total 

containing  information,  include 
Annexes,  Appendices,  etc.) 

41 

7.  DESCRIPTIVE  NOTES  (the  category  of  the  document,  e.g.  technical  report,  technical  note  or  memorandum.  If  appropriate,  enter  the  type  of 
report,  e.g.  interim,  progress,  summary,  annual  or  final.  Give  the  inclusive  dates  when  a  specific  reporting  period  is  covered.) 

DRDC  Ottawa  Technical  Memorandum 

8.  SPONSORING  ACTIVITY  (the  name  of  the  department  project  office  or  laboratory  sponsoring  the  research  and  development.  Include  the 
address.) 

Defence  R&D  Canada  -  Ottawa 
Ottawa,  Ontario,  Canada  K1 A  0Z4 


11.  DOCUMENT  AVAILABILITY  (any  limitations  on  further  dissemination  of  the  document,  other  than  those  imposed  by  security  classification) 

( X  )  Unlimited  distribution 

(  )  Distribution  limited  to  defence  departments  and  defence  contractors;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  defence  departments  and  Canadian  defence  contractors;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  government  departments  and  agencies;  further  distribution  only  as  approved 

(  )  Distribution  limited  to  defence  departments;  further  distribution  only  as  approved 

(  )  Other  (please  specify): 

12.  DOCUMENT  ANNOUNCEMENT  (any  limitation  to  the  bibliographic  announcement  of  this  document.  This  will  normally  correspond  to 
the  Document  Availability  (11).  However,  where  further  distribution  (beyond  the  audience  specified  in  11)  is  possible,  a  wider 
announcement  audience  may  be  selected.) 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


9a.  PROJECT  OR  GRANT  NO.  (if  appropriate,  the  applicable  research 
and  development  project  or  grant  number  under  which  the 
document  was  written.  Please  specify  whether  project  or  grant) 

15ec05 

10a.  ORIGINATOR'S  DOCUMENT  NUMBER  (the  official  document 
number  by  which  the  document  is  identified  by  the  originating 
activity.  This  number  must  be  unique  to  this  document.) 

DRDC  Ottawa  TM  2006-185 


9b.  CONTRACT  NO.  (if  appropriate,  the  applicable  number  under 
which  the  document  was  written) 


10b.  OTHER  DOCUMENT  NOS.  (Any  other  numbers  which  may 
be  assigned  this  document  either  by  the  originator  or  by  the 
sponsor) 


5.  DATE  OF  PUBLICATION  (month  and  year  of  publication  of 
document) 

September  2006 


6b.  NO.  OF  REFS  (total  cited  in 
document) 

19 


1.  ORIGINATOR  (the  name  and  address  of  the  organization  preparing  the  document. 
Organizations  for  whom  the  document  was  prepared,  e.g.  Establishment  sponsoring  a 
contractor's  report,  or  tasking  agency,  are  entered  in  section  8.) 

Defence  R&D  Canada  -  Ottawa 
Ottawa,  Ontario,  Canada  K1 A  0Z4 


2.  SECURITY  CLASSIFICATION 

(overall  security  classification  of  the  document, 
including  special  warning  terms  if  applicable) 

UNCLASSIFIED 


DCD03  2/06/87 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


13.  ABSTRACT  (a  brief  and  factual  summary  of  the  document.  It  may  also  appear  elsewhere  in  the  body  of  the  document  itself.  It  is  highly 
desirable  that  the  abstract  of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the 
security  classification  of  the  information  in  the  paragraph  (unless  the  document  itself  is  unclassified)  represented  as  (S),  (C),  or  (U). 

It  is  not  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is  bilingual). 


(U)  The  adaptive  local  polynomial  Fourier  transform  is  employed  for  the  improvement  of  the  ISAR  images  in  complex 
reflector  geometry  eases,  as  well  as  in  cases  of  fast  maneuvering  targets.  It  has  been  shown  that  this  simple  technique 
can  produce  significantly  improved  results  with  a  relatively  modest  calculation  burden.  Two  forms  of  the  adaptive 
LPFT  arc  proposed.  The  adaptive  parameter  in  the  first  form  is  calculated  for  each  radar  chirp.  An  additional 
refinement  is  performed  by  using  the  information  from  the  adjacent  chirps.  The  second  technique  is  based  on  the 
determination  of  the  adaptive  parameter  for  different  parts  of  the  radar  image.  The  numerical  analysis  demonstrates  the 
accuracy  of  the  proposed  techniques.  It  is  important  to  note  that  the  proposed  technique  docs  not  assume  any  particular 
model  of  radar  target  motion.  It  can  be  applied  for  any  realistic  motion  of  targets. 


14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a  document  and  could  be  helpful 
in  cataloguing  the  document.  They  should  be  selected  so  that  no  security  classification  is  required.  Identifiers  such  as  equipment  model 
designation,  trade  name,  military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords  should  be  selected  from  a 
published  thesaurus,  e.g.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  is  not  possible  to  select 
indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 

SAR 

ISAR 

Local  Polynomial  Fourier  Transform 

Time- Frequency  Analysis 

Moving  Targets 

Target  Detection 

Fourier  Transform 

Doppler  Processing 

Image  Analysis 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 

13.  ABSTRACT  (a  brief  and  factual  summary  of  the  document.  It  may  also  appear  elsewhere  in  the  body  of  the  document  itself.  It  is  highly 
desirable  that  the  abstract  of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the 
security  classification  of  the  information  in  the  paragraph  (unless  the  document  itself  is  unclassified)  represented  as  (S),  (C),  or  (U). 

It  is  not  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is  bilingual). 


(U)  The  adaptive  local  polynomial  Fourier  transform  is  employed  for  the  improvement  of  the  ISAR  images  in  complex 
reflector  geometry  cases,  as  well  as  in  cases  of  fast  maneuvering  targets.  It  has  been  shown  that  this  simple  technique 
can  produce  significantly  improved  results  with  a  relatively  modest  calculation  burden.  Two  forms  of  the  adaptive 
LPFT  are  proposed.  The  adaptive  parameter  in  the  first  form  is  calculated  for  each  radar  chirp.  An  additional 
refinement  is  performed  by  using  the  information  from  the  adjacent  chirps.  The  second  technique  is  based  on  the 
determination  of  the  adaptive  parameter  for  different  parts  of  the  radar  image.  The  numerical  analysis  demonstrates  the 
accuracy  of  the  proposed  techniques.  It  is  important  to  note  that  the  proposed  technique  does  not  assume  any  particular 
model  of  radar  target  motion.  It  can  be  applied  for  any  realistic  motion  of  targets. 


14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a  document  and  could  be  helpful 
in  cataloguing  the  document.  They  should  be  selected  so  that  no  security  classifcation  is  required.  Identifiers  such  as  equipment  model 
designation,  trade  name,  military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords  should  be  selected  from  a 
published  thesaurus,  e.g.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  is  not  possible  to  select 
indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 

SAR 

ISAR 

Local  Polynomial  Fourier  Transform 

Time-Frequency  Analysis 

Moving  Targets 

Target  Detection 

Fourier  Transform 

Doppler  Processing 

Image  Analysis 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


Defence  RfrD  Canada 


R  &  D  pour  la  defense  Canada 


Canada's  Leader  in  Defence  Chef  de  file  au  Canada  en  matiere 

and  National  Security  de  science  et  de  technologie  pour 

Science  and  Technology  la  defense  et  la  securite  nationale 


www.drdc-rddc.gc.ca 


